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ABSTRACT 

A new technique is presented for determining the black-hole masses of high-redshift 
quasars from optical spectroscopy. The new method utilizes the full-width half max- 
imum (FWHM) of the low-ionization Mgn emission line and the correlation between 
broad-line region (BLR) radius and continuum luminosity at 3000A. Using archival 
UV spectra it is found that the correlation between BLR radius and 3000A luminos- 
ity is tighter than the established correlation with 5100A luminosity. Furthermore, it 
is found that the correlation between BLR radius and 3000A continuum luminosity 
is consistent with a relation of the form Rblr oc Ai°' 5 , as expected for a constant 
ionization parameter. Using a sample of objects with broad-line radii determined from 
reverberation mapping it is shown that the FWHM of Mg n and H/3 are consis- 
tent with following an exact one-to-one relation, as expected if both H/3 and Mgn 
are emitted at the same radius from the central ionizing source. The resulting virial 
black-hole mass estimator based on rest-frame UV observables is shown to reproduce 
black-hole mass measurements based on reverberation mapping to within a factor of 
2.5 (lc). Finally, the new UV black-hole mass estimator is shown to produce identical 
results to the established optical (H/3) estimator when applied to 128 intermediate- 
redshift (0.3 < z < 0.9) quasars drawn from the Large Bright Quasar Survey and the 
radio-selected Molonglo quasar sample. We therefore conclude that the new UV virial 
black-hole mass estimator can be reliably used to estimate the black-hole masses of 
quasars from z ~ 0.25 through to the peak epoch of quasar activity at z ~ 2.5 via 
optical spectroscopy alone. 

Key words: galaxies: fundamental parameters - galaxies: active - galaxies: nuclei - 
galaxies: high redshift - quasars: general - quasars: emission lines 



1 INTRODUCTION 

A reliable method for estimating the black-hole masses of 
high-redshift quasars would provide crucial new information 
for understanding the nature and cosmological evolution of 
quasars. In this paper we address this problem by deriving 
a virial black-hole mass estimator based on rest-frame UV 
observables which is capable of producing reliable quasar 
black-hole mass estimates from optical spectra out to red- 
shifts of z ~ 2.5. 

Currently, the most direct measurements of the central 
black-hole masses of powerful active galactic nuclei (AGN) 
come from the reverberation mapping studies of 17 Seyferts 
and 17 PG quasars by Wandel, Peterson & Malkan (1999) 
and Kaspi et al. (2000) respectively. By measuring the time- 
lag between the variations of the AGN continuum and the 
broad H/3 (4861 A) emission line, these long-term monitor- 
ing programmes have delivered direct measurements of the 
radius of the H/3 emission region from the central ionizing 



source (Rblr = cr, where r is the time lag between the con- 
tinuum and H/3 variations). By combining these measure- 
ments with the assumption that the FWHM of the H/3 line 
reflects the virialized bulk motion of the line-emitting ma- 
terial, it is then possible to produce the so-called virial esti- 
mate of the central black hole mass (see Wandel, Peterson & 
Malkan 1999 for a full discussion) . Although the virial mass 
estimate is potentially subject to many uncertainties (eg. 
Krolik 2001), comparison with the predictions of the com- 
pletely independent relation between black-hole mass and 
host-galaxy stellar velocity dispersion have shown the two 
to be in excellent agreement (Ferrarese et al. 2001, Gebhardt 
et al. 2000). 

Unfortunately, because of the many years of monitor- 
ing required, it is unrealistic to expect that reverberation 
mapping measurements can be obtained for large samples of 
quasars, even at low redshift. However, the existing rever- 
beration mapping of low redshift AGN has revealed a corre- 
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lation between Rblr and the monochromatic AGN contin- 
uum luminosity at 5100A (eg. Rblr oc XL°i 00 ; Kaspi et al. 

2000) . By exploiting this correlation to estimate Rblr it is 
possible to produce a virial black-hole mass estimate from a 
single spectrum covering the H/3 emission line. 

This technique has recently been widely employed to 
investigate how the masses of quasar black holes relate to 
the properties of the surrounding host galaxies (eg. McLure 
& Dunlop 2002; 2001; Laor 2001) and the radio luminosity 
of the central engine (Dunlop et al. 2002; Lacy et al. 2001). 
However, because of the atmospheric restrictions and rela- 
tive observational expense of near-infrared spectroscopy, the 
vast majority of previous studies have relied on optical spec- 
troscopy of the H/3 region. Unfortunately, this effectively en- 
forces a redshift upper limit ofz~0.8 — 0.9, with the vast 
majority of studies concentrating on samples at z < 0.3. 

An ideal solution to this problem would be to cali- 
brate a rest-frame UV emission line and continuum mea- 
surement such that it becomes possible to produce a virial 
black-hole mass estimate for high-redshift quasars from rela- 
tively straightforward optical spectroscopy. It is the purpose 
of this paper to provide the required UV calibration. There 
are two strong, permitted, broad UV lines which are candi- 
dates to replace H/3 in the virial black-hole mass estimate; 
Civ(1549A) and Mg n(2798A). We note here that Vester- 
gaard (2002) recently derived a UV virial black-hole mass 
estimator using the Civ FWHM and the continuum lumi- 
nosity at 1350A . Although Civ does have the advantage of 
being accessible with optical spectroscopy at the very high- 
est redshifts (3 < z < 5), where Mgn is redshifted into the 
near-infrared, there are a number of reasons to believe that 
Mgn represents a better substitute for H/3 over the redshift 
range 0.8 < z < 2.5. 

The main reason for adopting Mgn as the UV tracer 
of BLR velocity is that, like H/3, Mgn is a low-ionization 
line. Furthermore, due to the similarity of their ionization 
potentials, it is reasonable to expect that the Mgn and H/3 
emission lines are produced by gas at virtually the same ra- 
dius from the central ionizing source, an assumption which 
is supported by the good agreement between the Mgn and 
H/3 FWHM (see Section 5). This has two important conse- 
quences. Firstly, it allows us to directly adopt the reverbera- 
tion mapping determinations of Rblr when calibrating the 
correlation between Rblr and 3000A continuum luminos- 
ity. Secondly, because the line- widths of Mgn and H/3 should 
trace the same BLR velocities, we are able to simply sub- 
stitute the FWHM of Mgn for that of H/3 in the virial mass 
estimator. 

There is one further practical advantage which favours 
Mgn over Civ as the UV substitute for H/3 . Due to the 
fact that Mgn becomes accessible to optical spectroscopy 
at z ~ 0.3, both H/3 and Mgn can be observed in a single 
optical spectrum within the redshift interval 0.3 < z < 0.9. 
Consequently, it is therefore possible to directly compare the 
results of the new UV virial mass estimator against those of 
the established optical virial mass estimator for objects at 
intermediate redshifts. This test is performed in Section ^ 
using data for the Large Bright Quasar Survey (Forster et al. 

2001) and the radio-selected Molonglo quasar sample (Baker 
et al. 1999). 

The structure of this paper can be summarized as fol- 
lows. In Section we briefly review the main aspects of 



the virial black-hole mass estimate, before proceeding in 
Section |^ to describe the properties of the sample of low- 
redshift AGN with reverberation mapping results. In Sec- 
tions ^ and ^ this sample is used to calibrate the relation 
between broad-line radius and 3000A continuum luminosity 
{Rblr — AL3000) and to demonstrate the viability of adopt- 
ing Mgn as a direct substitute for H/3 in the virial black-hole 
mass estimator. Our final calibration of the UV black-hole 
mass estimator is presented in Section ^. In Section [7] the 
effectiveness of the new UV black-hole mass estimator is 
demonstrated by application to members of the Large Bright 
Quasar Survey (LBQS) and the Molonglo quasar sample 
(MQS) for which both H/3 and Mgn FWHM measurements 
are available. Our conclusions are presented in Section ^[ 
All cosmological calculations presented in this paper assume 
n m = 0.3, A = 0.7, Ho = 70 kms _1 Mpc _1 . 



2 THE VIRIAL BLACK-HOLE MASS 
ESTIMATE 

The underlying assumption behind the method of estimat- 
ing AGN black-hole masses from the width of their broad 
emission lines is that the motion of the line-emitting ma- 
terial is virialized. Under this assumption the width of the 
broad lines can be used to trace the Keplerian velocity of the 
broad-line gas, and thereby allow an estimate of the central 
black-hole mass via the formula : 

M bh = G^RblrVIlr (1) 

where Rblr is the BLR radius and Vblr is the Keplerian 
velocity of the BLR gas. Recent evidence to support the 
Keplerian interpretation of AGN broad-line widths has been 
presented by Peterson & Wandel (2000) and Onken & Peter- 
son (2002). These authors demonstrate, at least for several 
Seyfert galaxies which have been well studied for reverbera- 
tion mapping purposes, that the FWHM of various emission 
lines follow the Vblr oc Rg LR relation expected for Keple- 
rian motion. It is important to note that it is not clear that 
the virial condition is satisfied for the high-ionization Civ 
emission line. Substantial evidence exists in the literature 
that Civ may be produced in some form of outflow, as is 
suggested by the systematic blue-shift of Civ compared to 
H/3 (eg. Marziani et al. 1996). 



2.1 Broad-line region geometry 

In Eqn |l] the velocity of the BLR gas is taken as 
Vblr = fx i?/3(FWHM), where / is a geometric factor 
which relates the H/3 FWHM to the intrinsic Keplerian ve- 
locity. Due to the fact that there is currently no consensus 
on the geometry of the BLR in radio-quiet quasars it is con- 
ventional to set / = v3/2. However, it is clear that if AGN 
broad lines are in general produced in a more disc-like con- 
figuration, evidence for which has been found in radio-loud 
quasars by numerous authors (eg. Wills & Browne 1986; 
Brotherton 1996; Vestergaard, Wilkes & Barthel 2000), then 
/ will be inclination dependent. Indeed, in McLure & Dun- 
lop (2002) evidence was presented that the H/3 FWHM is 
inclination dependent in radio-quiet quasars, and more con- 
sistent with the orbits of the H/3 emitting material having 
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a flattened disc-like geometry than being randomly orien- 
tated. However, even in the scenario where the FWHM of 
H/9 and Mgn are entirely orientation dependent, inclination 
will only have a significant effect on the virial black-hole 
mass estimate for those objects within ~ 15° of the line of 
sight. If quasars are taken to be distributed randomly within 
~ 45° of the line of sight (Barthel 1989) then this will only 
effect some 10% of objects. Consequently, for the purposes 
of providing a generally applicable UV black-hole mass es- 
timator we adopt / = 1 throughout this paper. 



3 THE REVERBERATION MAPPED SAMPLE 

As mentioned previously, the sample of 17 Seyferts and 
17 PG quasars with reverberation mapping measurements 
(hereafter the RM sample) currently represents the best 
available set of black-hole mass estimates of powerful broad- 
line AGN. Therefore, it is the objects from this sample which 
we will use to calibrate the new UV virial black-hole mass 
estimator. Optical and UV data for the combined 34-object 
RM sample are listed in Table [|. The Rblr data are taken 
from Kaspi et al. (2000), with the exception of NGC 4051 for 
which we take the updated figure of Rblr = 5.9 light-days 
from Peterson et al. (2000). 

3.1 Optical data 

The 5100A luminosities for the PG quasars have been calcu- 
lated using the fluxes from the spectrophotometric study of 
Neugebauer et al. (1987). For the Seyfert galaxies the 5100A 
luminosities are from Kaspi et al. (2000) after conversion to 
our cosmology. The H/3 FWHM values for the PG quasars 
are taken from Boroson & Green (1992). The H/3 FWHM 
values for the Seyfert galaxies are drawn from a number of 
literature sources which are listed in the caption to Table [j] 

3.2 UV data 

As with the optical continuum luminosities, the 3000A lu- 
minosities listed in Table |l| for the PG quasars are derived 
from Neugebauer et al. (1987). In contrast, the 3000A lumi- 
nosities for the Seyfert galaxies, together with 17/22 of the 
Mgn FWHM values, are new measurements based on our fit- 
ting of spectra from the International Ultraviolet Explorer 
(IUE) final data archive. A full discussion of our line-fitting 
method will be presented in a follow-up paper (McLure et 
al. 2002, in prep). In brief, we fit the region spanning the 
Mgn emission line (2300A^ 3100A) with a combination of 
a power-law continuum, a Fell emission template and two 
gaussians representing the broad and narrow components of 
the Mgn emission line. 

As with the H/3 emission line, it is crucial to account 
for blending with broad Fen emission in order to obtain a 
reliable measurement of the Mgn line-width. Consequently, 
following the analysis of Corbin & Boroson (1996) we have 
constructed a template of the Fell emission around the Mgn 
line from an archival IUE spectrum of the narrow line Seyfert 
galaxy I Zw 1. During the fitting process the amplitude and 
FWHM of the Fell template were left as free parameters in 
order to obtain the best possible match to the Fell emission 
of each object. 



4 ESTIMATING THE BLR RADIUS 

In order to provide an estimate of Rblr it is necessary, in 
the absence of reverberation mapping, to exploit the known 
correlation between Rblr and AGN continuum luminosity. 
By combining the results of their reverberation mapping ob- 
servations with those of Wandel, Peterson & Malkan (1999), 
Kaspi et al. (2000) investigated the correlation between 
Rblr and 5100A continuum luminosity and found a best-fit 
of the form Rblr oc AL^qq 03 . In their re-analysis of the 
same correlation, Peterson et al. (2000) found a best fitting 
relation of the form Rblr oc ALgjQo ±0 ' 02 . As well as being 
important within the context of the virial black-hole mass 
estimate, the slope of the Rblr — ALa relation is of intrin- 
sic interest because it provides information about the AGN 
ionization parameter U (U = j-^r , where Q = J j^jdv). 
For example, neither of the two fits mentioned above are 
consistent with the relation Rblr oc AL^' 5 expected if the 
ionization parameter is approximately the same for all AGN, 
regardless of luminosity. In this section we investigate the 
slope of the Rblr — \L\ relation using continuum luminos- 
ity measurements at both 3000A and 5100A. 

4.1 Regression analysis 

Given the potential influence of the Rblr — AL^ relation 
upon the accuracy of the virial black-hole mass estimator, 
we have chosen to employ three different regression tech- 
niques in our analysis to explore the possible variation in 
the best-fitting slope. The first technique is a straightfor- 
ward weighted least-squares fit (WLS), in which the fit is 
weighted by the error in the reverberation mapping radius 
(Press et al. 1992). The second technique is a robust itera- 
tive chi-square fit (FITxy) in which the errors in both vari- 
ables are taken into account (Press et al. 1992). However, 
although FITxy takes into account the errors on both vari- 
ables, it does not account for the fact that there is likely to 
be intrinsic scatter in the Rblr — A£a relation. To account 
for this, the third technique we have employed is the BCES 
estimator of Akritas & Bershady (1996) which accounts for 
errors in both variables and intrinsic scatter. 

4.2 The radius-luminosity relation at 3000A and 
5100A 

The results of the fitting process (in log — log space) of both 
the Rblr — AL3000 and Rblr — AL5100 relations are listed 
in Table ^ where the listed BCES fit is the bisector. It can 
be seen from Table ^ that the results of all three regression 
methods are consistent with each other for both the Rblr — 
AL5100 and Rblr — AL3000 relations. Furthermore, it can 
be seen that all three regression methods indicate that the 
Rblr — AL3000 relation is fully consistent with a slope of 
AL°' 5 , as expected for a constant ionization parameter. 

The three fits to the Rblr — AL5100 relation indicate 
that a slope of a AL°' 6 is more appropriate for our choice of 
continuum data, in good agreement with the fit of Peterson 
et al. (2000). This is flatter than the slope of AL°' 7 previously 
determined for the Rblr ~ AL5100 relation by Kaspi et al. 
(2000), although the two are formally consistent. 

There are two effects which are responsible for our 
finding of a flatter slope in the Rblr — AL5100 relation. 
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Object 


z 


log (AL3000) 
/W 


log (AL5100) 
/W 


log (Rblr) 
/lt-days 


Mgn FWHM 
/kms -1 




H/3 FWHM 
/kms -1 




PG 0026+129 


0.142 


38.23 


38.07 


2.05 






1860 


1 


PG 0052+251 


0.155 


38.29 


38.04 


2.13 






5200 


1 


PG 0804+761 


0.100 


38.40 


37.96 


2.19 






3070 


1 


PG 0844+349 


0.064 


37.61 


37.44 


1.38 


3800 


2 


3800 


1 


PG 0953+414 


0.239 


38.69 


38.41 


2.18 






3130 


1 


PG 1211+143 


0.085 


38.20 


38.00 


2.00 


2050 


5 


1860 


1 


PG 1226+023 


0.158 


39.23 


38.97 


2.59 


3160 


5 


3520 


1 


PG 1229+204 


0.064 


37.53 


37.33 


1.70 


2560 


5 


3360 


1 


PG 1307+085 


0.155 


38.26 


38.01 


2.09 






5320 


1 


PG 1351+640 


0.087 


37.93 


37.79 


2.36 


2790 


5 


5660 


1 


PG 1411+442 


0.089 


37.68 


37.53 


2.01 


1940 


5 


2670 


1 


PG 1426+015 


0.086 


38.11 


37.78 


1.98 






6820 


1 


PG 1613+658 


0.129 


37.85 


37.74 


1.59 






8450 


1 


PG 1617+175 


0.114 


38.16 


37.74 


1.93 






5330 


1 


PG 1700+518 


0.292 


38.75 


38.70 


1.94 






2210 


1 


PG 1704+608 


0.371 


38.84 


38.72 


2.50 






6560 


1 


PG 2130+099 


0.061 


37.77 


37.47 


2.30 


4360 


2 


2330 


1 


3C 120 


0.035 


37.80 


36.94 


1.62 


3360 


2 


1910 


3 


3C 390 


0.052 


36.49 


36.89 


1.36 


8410 


2 


10000 


3 


Akn 120 


0.035 


37.60 


37.22 


1.57 


9660 


2 


5800 


3 


F9 


0.047 


37.62 


37.22 


1.21 


5680 


2 


6080 


4 


IC 4329A 


0.016 


34.25 


36.28 


0.15 






5050 


3 


Mrk 79 


0.022 


36.42 


36.70 


1.26 


4730 


2 


4470 


3 


Mrk 110 


0.035 


36.77 


36.66 


1.27 


2420 


2 


2120 


1 


Mrk 335 


0.026 


37.15 


36.87 


1.21 


8605 


2 


5685 


6 


Mrk 509 


0.033 


37.67 


37.24 


1.88 


4790 


2 


2270 


1 


Mrk 590 


0.026 


36.88 


36.78 


1.30 


4370 


2 


2470 


3 


Mrk 817 


0.032 


36.99 


36.80 


1.19 


5770 


2 


4490 


3 


NGC 3783 


0.009 


36.32 


36.31 


0.65 


2690 


2 


3780 


4 


NGC 3227 


0.004 


35.18 


35.37 


1.04 






4920 


3 


NGC 4051 


0.003 


35.08 


34.78 


0.77 


2790 


2 


1170 


3 


NGC 4151 


0.003 


34.93 


35.92 


0.48 


3990 


2 


5910 


3 


NGC 5548 


0.014 


36.52 


36.50 


1.33 


5420 


2 


6200 


4 


NGC 7469 


0.017 


36.95 


36.81 


0.69 


3258 


2 


3000 


3 



Table 1. Data for the reverberation mapped sample. Columns 1 and 2 detail the name and redshift of each object. Columns 3 and 4 list 
the monochromatic luminosities at 3000A and 5100A respectively. The 3000A and 5100A luminosities for the PG quasars are derived 
from the spectrophotometric study of Neugebauer et al. (1987). For the Seyfert galaxies the 3000A luminosities are derived from new 
fits to IUE archive spectra, while the 5100A luminosities are taken from Kaspi et al. (2000) after conversion to our cosmology. Column 5 
lists the BLR radius measurements from Kaspi et al. (2000), expect for NGC 4051 which is taken from Peterson et al. (2000). Columns 
6 & 7 detail the FWHM of the Mgn and H/3 emission lines respectively. The references for the FWHM measurements are as follows : 1. 
Boroson & Green (1992), 2. This work, 3. Peterson, Wandel & Malkan (1999), 4. Marziani et al. (1996), 5. Corbin & Boroson (1996) and 
6. Zheng et al. (1995) 



Firstly, our choice of cosmology (Q m = 0.3, A = 0.7, Ho = 
70 kms -1 Mpc -1 ) results in higher luminosities for the 17 
PG quasars in the RM sample, compared to the Q m = 1, 
A = 0, Ho — 75 kms _1 Mpc _1 cosmology adopted by Kaspi 
et al., while having little effect on the luminosities of the 
lower redshift Seyfert galaxies. Secondly, further flattening 
of the slope is due to the fact that the 5100A fluxes for the 
PG quasars adopted here, derived from the Neugebauer et 
al. (1987) data, are typically ~ 0.1 dex higher than those 
determined by Kaspi et al. Indeed, if we adopt the same PG 
quasar luminosities and cosmology as Kaspi et al., then the 
BCES bisector fit to the Rblr — AL5100 relation becomes: 
Rblr oc ALgjoo 3 entirely consistent with the slope de- 
termined by Kaspi et al. 

Throughout the remainder of the paper we adopt the 
BCES bisector fits as our best estimate of the Rblr — ALa 
relation at both 3000 A and 5100A . The bisector fits to both 



data-sets are shown in Fig [l| and have linear versions as 
follows: 

Rblr = (25.2 ± 3.0) [AL 3 ooo/10 37 W] (0 - 47±a05) (2) 

Rblr = (26.4 ± 4.4) [AL 5100 /10 37 ( °' 61±ai0) (3) 
where Rblr is in units of light-days. 



5 ESTIMATING THE BROAD-LINE 
VELOCITY 

After successfully calibrating the Rblr — AL3000 relation, 
the second required element of the UV black-hole mass esti- 
mator is to calibrate the Mgn FWHM as a substitute for 
H/3 . As discussed in Section [j], the principle reason for 
adopting Mgn as the UV H/3 proxy is that it is a strong, 
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Log(AL x (3000)/W) Log(XL x (5100)/W) 

Figure 1. The correlation between broad-line radius and AGN continuum luminosity at 3000A and 5100A . The 17 PG quasars are 
shown as filled squares and the 17 Seyfert galaxies are shown as filled circles. The best-fitting BCES bisector fits to both correlations are 
shown as solid lines, and correspond to Rblr ^3000 an< ^ Rblr oc AL^qq respectively. 



Method 


A 


b 


a 


WLS 


5100A 


0.57 ±0.02 


-19.53 ±0.80 


FITxy 


5100A 


0.60 ± 0.02 


-20.75 ±0.92 


BCES 


5100A 


0.61 ± 0.10 


-21.00 ±3.67 


WLS 

FITxy 

BCES 


3000A 
3000A 
3000A 


0.50 ±0.02 
0.53 ± 0.02 
0.47 ± 0.05 


-17.19 ±0.70 
-18.08 ±0.82 
-16.10 ± 1.69 



Table 2. Results of the regression analysis of the correlation be- 
tween BLR radius and the AGN continuum luminosity at 5100A 
and 3000A . The fits were performed in log — log space and are of 
the form: log Rblr = blog(\L\) + a. The three different regres- 
sion methods are discussed in the text. 



fully-permitted, low-ionization line, and as such should pro- 
vide the best UV analog of H/3 . Moreover, because the ion- 
ization potentials of H/3 and Mgn are very similar, the two 
lines should both be emitted at approximately the same ra- 
dius from the central ionizing source. Fortunately, this as- 
sumption can be tested directly using the 22 objects from the 
RM sample for which we have Mgn FWHM measurements. 
If Mgn and H/3 are being produced at the same broad-line 
radius then, if the BLR is virialized, we should see a 1:1 re- 
lation between Mgn FWHM and H/3 FWHM. Fig g shows 
Mgn FWHM versus H/3 FWHM for the 22 objects from 
the RM sample with Mgn FWHM measurements. The solid 
line in Fig H shows an exact 1:1 relation between Mgn and 
H/3 FWHM. It can be seen that, with the exception of NGC 
4051, the objects in the RM sample are perfectly consis- 
tent with Mgn and H/3 tracing the same BLR velocity. The 
clear outlier, NGC 4051, is a well studied example of a 
narrow-line Seyfert 1 (NLS1). The H/3 FWHM for this ob- 
ject has been determined by numerous authors and is consis- 
tently found to be around 1000 kms -1 (eg. Wandel, Peter- 
son & Malkan (1999) find H/3 FWHM=1170 kms" 1 ). How- 
ever, from our line-fitting of IUE archive spectra we find 



a strong broad component to the Mgn emission line with 
FWHM=2790 kms" 1 . 

In fact, several previous studies have found that the UV 
emission lines of NLS1 possess strong broad components 
which are not seen in their Balmer lines (eg. Rodrfguez- 
Pascual, Mas- Hesse & Santos-Lleo 1997, Zheng et al. 1995). 
This suggests that Mgn may actually provide a better es- 
timate of the BLR velocity of low-ionization gas in NLS1, 
where the Balmer lines could be biased, at least in part, 
by inclination effects. Irrespective of the reason why NGC 
4051 is an outlier in Fig ^, we can safely exclude it from 
a regression fit to the Mgn - H/3 FWHM relation. The rea- 
son for this is that we are principally interested in providing 
a UV black-hole mass estimator for high redshift quasars, 
the vast majority of which have H/3 FWHM > 2000 kms -1 . 
The dotted line in Fig g shows the BCES bisector fit which 
is equivalent to a relation of the form: 

H/3(FWHM) oc Mg//(FWHM) 102±0 - 14 (4) 

which is clearly consistent with a linear relation. In the light 
of this result, an exact 1:1 relation between the Mgn FWHM 
and H/3 FWHM is assumed in the next section for the final 
calibration of the UV black-hole mass estimator. 



6 THE UV BLACK-HOLE MASS ESTIMATOR 

Having determined in the previous two sections the Rblr — 
AI/3000 relation and the 1:1 scaling between H/3 FWHM and 
Mgn FWHM, we are now in a position to derive our UV 
black-hole mass estimator. In Fig |§] we show the reverbera- 
tion mapping black-hole mass estimate versus the new UV 
mass estimate, which uses the Rblr — AL3000 relation to 
estimate Rblr and the Mgn FWHM to trace the BLR ve- 
locity. It is clear from Fig that Mrk 335 is an outlier. As 
with NGC 4051, Mrk 335 is another example of a NLS1. The 
reason Mrk 335 is an outlier in Fig | is that the FWHM 
of the variable component of its H/3 line determined from 
reverberation mapping is only 1260 kms -1 (Wandel, Pe- 
terson & Malkan 1999), while from our line- fitting of IUE 
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Figure 2. A Log-Log plot of Mgn FWHM versus H/3 FWHM 
for the 22 objects from the RM sample for which it was possible 
to obtain measurements of both line- widths. The solid line is an 
exact 1:1 relation. The dotted line is the BCES bisector fit, ex- 
cluding NGC 4051, and has slope of 1.02 ± 0.14. A 15% error in 
the FWHM of H/3 and Mgn has been assumed. Symbols are as 
Fig|. 

archive spectra we detect a strong broad component in the 
Mgn line with FWHM=8605 kms -1 . Our line-fitting results 
for this object are in excellent agreement with those of 
Zheng et al. (1995) who detected a broad Mgn component 
of FWHM=8218 kms -1 . Again, as mentioned in the previ- 
ous section, it appears a possibility that for at least some 
NLSl's the width of the H/3 emission line may not be a re- 
liable tracer of BLR velocities. Of course, it is also possible 
that Mrk 335 may simply be a genuine outlier. The solid 
line in Fig ^ is the bisector fit to the data, excluding Mrk 
335, and has the form: 



log M bh (RM) = 1.12(±0.22) \ogM bh {UV) 
-1.08(±1.75) 



(5) 



again, consistent with a linear relation. Consequently, our 
final calibration of the UV mass estimator has the form: 



log M bh (RM) = log M bh (UV) - 0.16 



(6) 



where the off-set of 0.16 is adopted to ensure a mean 
M b h{RM):M b h(UV) ratio of unity. In terms of a useful for- 
mula the best fitting calibration of the UV black-hole mass 
estimator is therefore: 



Mbh 
M 



3.37 



/ AI/300Q \ 1 

U0 37 W/ 



/ FWHM(MgII) Y 
\ kms" 1 J 



(T) 



For completeness we have also re-derived the optical (H/3) 
virial black-hole mass estimator using our re-analysis of the 
Rblr — AL5100 relation in Section Id As expected, the cor- 
relation between M b h{RM) and M b h(Hf3) is also consistent 
with being linear (M bh (RM) oc M bh (H /3)°- 87±0 - 18 ) which, 
in the same fashion as for the UV estimate, leads to: 



M b h _ ^ 74 / ALsiop y 
VI'.; " 



V10 37 W7 



/ FWHM(Hf3) X 
\ kms" 1 ) 



(8) 




Log (M bh /M )(UV) 

Figure 3. Full reverberation mapping black-hole mass estimate 
(based on Rblr measurements and the rms H/3 FWHM) plotted 
against the UV mass estimate (based on the Rblr — ^3000 
relation of Eqn 2 and the FWHM of Mgn) for the 22 objects 
from the RM sample for which it was possible to obtain Mgn 
FWHM measurements. The solid line is the BCES bisector fit, 
excluding Mrk 335, which has a slope of 1.12 ± 0.22. The dotted 
line is the adopted linear relation (Eqn 6). Symbols as in Fig |l|. 



6.1 Accuracy of the UV estimator 

Having arrived at the final calibration of the UV black-hole 
mass estimator it is of obvious interest to assess the accu- 
racy with which it can reproduce the full reverberation mass 
estimate. Excluding Mrk 335, the mean difference between 
the reverberation and the UV estimator is : 



< \ag(M bh )(RM) - \og(M bh )(UV) >= 0.00 ± 0.40 



(9) 



where the quoted uncertainty is the standard deviation (a) 
and not the standard error. Therefore, provided the RM 
sample is representative of broad-line AGN in general, we 
conclude that the UV black-hole estimator provided by 
Eqn can reproduce the reverberation black-hole mass to 
within a factor of 2.5 (lcr). For the same sample the uncer- 
tainty on the optical estimator given in Eqn H is a factor of 
2.7. 



7 APPLICATION TO QUASAR SURVEYS: 
THE LBQS AND MQS 

Having derived the desired black-hole mass estimator in the 
previous section, we next assess its performance when ap- 
plied to the spectra of the LBQS and the radio-selected 
MQS. If the new UV black-hole mass estimator is to be 
used in the analysis of high-redshift quasars, it is clearly 
important to ensure that it does not produce results which 
are biased with respect to the usual optical (H/3) estima- 
tor. Therefore, our intention in this section is to simply test 
how the new UV mass estimator (Eqn [?]) compares to the 
usual optical mass estimator (Eqn |8|) when applied to those 
objects from the LBQS and MQS for which both the ap- 
propriate UV (AL3000, Mgn FWHM) and optical (AL5100, 
H/3 FWHM) data are available. 
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Due to the fact that the LBQS is optically selected, 
while the MQS is a complete radio-selected quasar sample, 
applying the UV black-hole mass estimator to both samples 
allows us to ensure that Eqn [j] is equally applicable to both 
radio-quiet and radio-loud quasars. 



7.1 The Large Bright Quasar Survey 

The LBQS consists of 1058 optically-selected quasars in the 
redshift range 0.2 < z < 3.3 (see Hewett, Foltz & Chaf- 
fee (1995) for a full description). A comprehensive study 
of the optical continuum and emission line properties of 
992 quasars from the LBQS has recently been published by 
Forster et al. (2001). In this section we use data from Forster 
et al. for a 99-object sub-sample with reliable measure- 
ments of both Mgn and H/3 FWHM and continuum fluxes 
at 3000A and 5100A . This sub-sample comprises 68% of 
the 145 LBQS objects in the redshift range 0.20 < z < 0.66 
within which both Mgn and H/3 fall on the LBQS spectra. 



7.2 The Molonglo quasar sample 

The Molonglo quasar sample is a complete, radio-selected, 
sample consisting of 111 quasars in the redshift range 
0.1 < z < 2.9 with a 408 MHz flux density greater than 
0.95 Jy (Kapahi et al. 1998). Optical spectra for 79 MQS 
quasars were published by Baker et al. (1999), together with 
a preliminary analysis of their emission line properties. We 
are currently engaged in an automated study of the optical 
and UV emission lines from the MQS optical spectra (kindly 
provided by J. Baker), from which it has been possible to 
reliably measure the Mgn and H/3 FWHM for 29/38 objects 
where the MQS spectra include both emission lines. Unlike 
the LBQS, Baker (1997) has shown that many MQS ob- 
jects display substantial reddening in their optical spectra. 
As a result, the continuum luminosities of the reddest MQS 
objects were de-reddened using the SMC extinction curve of 
Pei (1992), assuming an intrinsic continuum slope of a — 0.5 
(fv oc v~ a ) which was found by Baker (1997) to be typical 
of the least reddened MQS objects. 



7.3 Optical versus UV black- hole mass estimators 

In Fig ^ we show the optical black-hole mass estimator plot- 
ted against the new UV black-hole mass estimator for a com- 
bined sample of 150 objects, comprising 99 from the LBQS, 
29 from the MQS and the 22 objects of the RM sample. 
Also shown is the BCES bisector fit to the LBQS and MQS 
objects which has the form: 



log M bh (Hf5) = 1.00(±0.08) log M bh (MgII) 
+0.06(±0.67) 



(10) 



which, as expected, is perfectly consistent with a linear re- 
lation. In Fig ^ we show a histogram of log M b u (Mgll) — 
log M bh (Hf3) for the 128 objects from the LBQS and MQS. 
The solid line shows the best-fitting gaussian which has 
a — 0.41. These results, in combination with those of Sec- 
tion ^, lead us to conclude that compared to the traditional 
optical black-hole mass estimator, the new UV estimator 
provides results which are unbiased and of equal accuracy. 




10 



Log(M b ./M ) Mgll 



Figure 4. The optical (H/3) versus UV (Mgn) virial black-hole 
estimators for 150 objects from the RM (filled circles), LBQS 
(filled squares) and MQS (open circles) samples. The solid line is 
the BCES bisector fit to the 128 objects from the MQS and LBQS 
samples and has a slope of 1.00 ±0.08. The outlier NGC4051 has 
been highlighted (see Section 5). 




Log(M bh (MgII)/M bh (H|S)) 

Figure 5. Histogram of log M bh (MgII) - log M bf JHf3) for the 
128 objects from the LBQS and MQS shown in Fig H. Also shown 
is the best-fitting gaussian which has a = 0.41. 



7.4 Black-hole mass and the radio-loudness 
dichotomy 

Although a full analysis of the black-hole masses of the 
LBQS and MQS samples is beyond the scope of this paper, 
we will briefly comment here on one obvious feature of Fig 4. 
It is immediately apparent from Fig ^ that the MQS quasars 
are confined to the highest black-hole masses. In fact, 27/29 
of the MQS quasars have black-hole masses > 1O 8 ' 5 M0, 
adopting either the UV or optical mass estimator. In con- 
trast, only 17% of the 99 LBQS quasars have UV and op- 
tical black-hole mass estimates which are both > 10 8 ' 5 M©. 
Naively this would appear to suggest that there is a clear 
division in black-hole mass between the radio-selected MQS 
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and the largely radio-quiet LBQS. However, because the 
wavelength coverage of the LBQS spectra only extends to 
~ 7500A, compared to ~ 10000A for those of the MQS, 
the need to include both Mgn and H/3 on the spectra intro- 
duces a considerable redshift bias into the two sub-samples 
under investigation here. Due to the differences in wave- 
length coverage the mean redshift of the LBQS sub-sample 
is 0.33 ± 0.01, while that of the MQS is 0.63 ± 0.03. Conse- 
quently, the members of the MQS sub-sample are on average 
more optically luminous than those of the LBQS sub-sample, 
with an average luminosity of log(AI/5ioo/W) = 38.14±0.08 
compared to log(AL 5 i o/W) = 37.64±0.07. Therefore, with 
reference to the Rblr — ALa relation, we would expect that, 
simply due to this luminosity bias, the average black-hole 
mass of the MQS sub-sample would be a factor of ~ 2 larger 
than that of the LBQS. However, in reality the average 
black-hole mass of the MQS sub-sample is 10 910±o or M 
whereas that of the LBQS sub-sample is more than a factor 
of 6 lower at l 8 27±0 06 M fl 

The origin of the additional factor of three difference in 
the average black-hole masses is a highly significant differ- 
ence in the distributions of H/3 and Mgn FWHM. For exam- 
ple, the average Mgn FWHM of the MQS sub-sample is a 
factor of 1.56±0.11 greater than the average Mgn FWHM of 
the LBQS sub-sample. The equivalent factor for H/3 FWHM 
is 1.68 ± 0.14. The application of a Kolmogorov-Smirnov 
test returns a probability of p = 1.33 x 10 -6 that the two 
H/3 FWHM distributions are drawn from the same parent 
distribution (p = 4.1 x 10~ 9 for the Mgn FWHM distribu- 
tions). 

It is important to note that the differences in the 
FWHM distributions of the two samples do not appear to 
be related to the previously mentioned redshift /luminosity 
bias. An application of the partial spearman rank correla- 
tion test (Macklin 1982) reveals no significant correlations 
between H/3 FWHM, redshift and AL5100 or Mgn FWHM, 
redshift and XLsooo- In fact, these results are in good agree- 
ment with previous studies of the black-hole masses of radio- 
quiet and radio-loud quasars which were not subject to 
a redshift /luminosity bias. In their study of a sample of 
radio-loud and radio-quiet quasars with a matched redshift- 
luminosity distribution at 2 < 0.5, McLure & Dunlop (2002) 
found the average black-hole mass of the radio-loud quasars 
to be a factor of ~ 2 larger than their radio-quiet coun- 
terparts, albeit with a large overlap. Likewise, Laor (2000) 
found a clean separation in black-hole mass within the PG 
quasar sample, with the radio-loud quasars being confined 
to Mbh ^ 10 9 Mq. In contrast, from their study of the First 
Bright Quasar Survey (FBQS), Lacy et al. (2001) concluded 
that the apparent gap in the quasar black-hole mass - radio- 
power plane is in fact filled by previously undetected radio- 
intermediate quasars. 

In either case, it is clear that genuinely powerful, 
low-frequency selected, radio-loud quasars harbour central 
black-hole masses drawn from the extreme end of the AGN 
black-hole mass function. At present it is not clear whether 
the low frequency of radio-quiet quasars with Mbh > 



* These figures are for the optical (H/3) mass estimator although, 
as expected from Fig ^, the black-hole mass difference is compa- 
rable if the UV mass estimator is adopted (AMn = 0.71 dex) 



10 8 ' 5 Mq is real or, alternatively, is due to an inclination 
selection effect whereby the line-widths of optically selected 
quasars are biased to low values by a preference for select- 
ing objects close to the line of sight. A more detailed analy- 
sis of the black-hole masses of high redshift radio-loud and 
radio-quiet quasars, the correlation between black-hole mass 
and radio power and the evidence for the role of inclination 
effects will be presented in a series of forthcoming papers 
(McLure et al. 2002, in prep). 

Finally, we note that our finding that the FWHM of 
H/3 and Mgn in the powerful radio-selected MQS objects 
are exclusively restricted to > 4000 kms -1 is of interest 
in the context of recent work on the location of radio- 
loud and radio-quiet AGN in so-called eigenvector 1 space 
(Marziani et al. 2001). In their study Marziani et al. con- 
cluded that powerful radio-loud quasars appear to be re- 
stricted to H/3 FWHM > 4000 kms" 1 , in excellent agree- 
ment with our analysis of the MQS spectra. 



8 CONCLUSIONS 

A new technique for estimating the central black- hole masses 
of high-redshift quasars using the Mgn FWHM and 3000A 
quasar continuum luminosity has been presented. The new 
technique has been calibrated using a sample of 34 low red- 
shift AGN with black-hole mass measurements based on 
long-term reverberation mapping experiments. The reliabil- 
ity of the new technique, with respect to the established 
optical virial black-hole mass estimator, has been tested us- 
ing published data for the LBQS together with the results 
of a new analysis of the emission line properties of the MQS. 
The main conclusions of this study can be summarized as 
follows: 

• The correlation between Rblr and monochromatic 
3000A continuum luminosity is found to display less scatter 
than that between Rblr and 5100A monochromatic con- 
tinuum luminosity. 

• The correlation between Rblr and 3000A continuum 
luminosity is found to be consistent with a relation of the 
form Rblr oc AL 9,5 , as expected for a constant ionization 
parameter in all AGN, irrespective of luminosity 

• The FWHM of Mgn is found to be an effective substi- 
tute for the FWHM of H/3 . The relationship between the 
two FWHM is found to be perfectly consistent with an exact 
1:1 scaling. 

• Combining the Rb lr~ AL3000 relation with the FWHM 
of Mgn produces a virial black-hole mass estimator based on 
rest-frame UV observables which is capable of reproducing 
black-hole masses determined from reverberation mapping 
to within a factor of 2.5 (la) 

• An application to objects from the LBQS and MQS 
demonstrates that the new UV black-hole mass estimator 
produces results which are unbiased, and of equal accuracy 
to the established optical (H/3) black-hole mass estimator. 

• We therefore conclude that the new UV black-hole mass 
estimator is ideal for determining quasar black-hole masses 
in the redshift range 0.25 < z < 2.5 via optical spectroscopy 
alone. 
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